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STOCHASTIC ANALYTICAL SOLUTION TO QUANTIFY THE EARTH'S 
SUBSURFACE MEAN HEAT FLOW AND ITS ERROR BOUNDS 



0 Field of the invention 

5 The present invention relates to obtaining the analytical closed form solution to mean 

and variance in heat flow by solving the stochastic heat conduction equation incorporating 
randomness in the thermal conductivity The method of the present invention has a wide 
range of applications in quantifying the thermal state of the crust and in obtaining closed 
form expressions for subsurface heat flow structure along with its error bounds. The exact 

10 formulae used in the invention can be used to better evaluate the thermal state for related oil 
and natural gas applications and also in tectonic studies and in studies related to the 
crystallization of minerals 
Background of the invention 

Heat flow studies play an important role in understanding the state of the lithosphere. 

1 S For modeling the conductive heat transfer in the crust the important controlling parameters 
are the radiogenic heat production and the thermal conductivity. In case of a stabilized 
continental crust the conductive heat transfer in steady state condition is a reasonably good 
approximation. Several authors have studied the evolution of the basal heat flow from a 
deterministic point of view. However due to the heterogeneous nature of the crust researchers 

20 are now solving the heat conduction problem from a stochastic point of view 

Quantification of perturbations in the temperatures and heat flow using stochastic 
analytical and random simulation techniques have been carried out by several authors. The 
uncertainty in the heat flow using a least squares inversion technique incorporating 
uncertainties in the temperature and thermal conductivities has been done, Tectonophysics, 

25 Vol 121, 1985 by Vausser et al . The effect of variation in heat source on the surface heat 
flow has also been studied, Journal Geophysical Research, V 91, 1986, by Vasseur and 
Singh, Geophysical Research Letters, VI 4, 1987, by Nielsen 

The small perturbation method to solve the stochastic heat conduction equation has 
been used to solve the l-D steady state equation with uncertainties in the heat sources to 

30 obtain the mean temperature field along with its error bounds , Geophysical Journal 
International, 135, 1998, by Srivastava and Singh. Several researchers have been using the 
numerical method, the random simulation method, to model the error structure in the thermal 
field incorporating uncertainties in the controlling thermal parameters, Tectooophysics, 
VI 56, 1988 by Royer and Danis, Marine and Petroleum Geology, V 14, 1997, by Gallagher 

35 et al, Tectonophysics, V 306, 1999a,b, by Jokinen and Kukkonen 
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The diffusion problems in stochastic framework are now being solved by yet another 
approach called the decomposition method, Journal of Hydrology, V 169, 1995, by Serrano 
In this invention the same approach as given in Geophysical Journal International, V 138, 
1 999, by Srivastava and Singh, has been used to solve the stochastic heat conduction equation 
5 incorporating Gaussian uncertainties in the thermal conductivity to obtain the solution to the 
mean and variance in the subsurface heat flow field. 
Objects of the invention 

The main object of the invention is to provide a method for the stochastic analysis to 
quantify the earth's subsurface area heat flow and its error bounds and which obviates the 
10 drawbacks of the prior art enumerated above. 

It is another object of the invention to provide a method for obtaining a solution to the 
stochastic heat conduction equation in an efficient manner which also enables the 
quantification of earth's subsurface area heat flow and its error bounds in an efficient manner. 
It is a further object of the invention to provide a method for the quantification of the 
15 thermal state of the earth* s crust wherein the variation in the thermal conductivity is 
quantified using a random correlation structure. 
Summary of the invention 

The present invention accordingly provides a method for the quantification of the 
earth's subsurface area heat flow and its error bounds by stochastic analysis comprising 
20 devising a stochastic heat conduction equation based on random thermal conductivity, an 
exponentially decreasing heat source and associated boundary conditions and arriving at a 
stochastic solution to the temperature field obtained using a series expansion method, and 
obtaining the expression for mean heat flow and variance m heat flow. 

In one embodiment of the invention, the boundary conditions are surface temperature 
25 and surface heat flow 

In another embodiment of the invention, the stochastic heat conduction equation is of 
the formula 

— K(z) — \ = -A(z) 
dz\ dz J 

In another embodiment of the invention, the stochastic heat conduction equation is 
30 solved using a series expansion method to obtain the closed form solution to the mean and 
variance in the heat flow fields, and using thermal conductivity of the subsurface layer of the 
earth's crust as a random parameter. 
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In yet another embodiment of the invention, the expression for heat flow is obtained 
and the expressions for the mean and variance in heat flow derived by taking expectation and 
using the property of the random field. 
Brief description of the accompanying drawings 

Figure I shows the plot of mean heat flow ± l. standard deviation versus depth for a 
synthetic model. 

Figure 2 shows the plot of mean heat flow ± 1 standard deviation versus depth for 
Latur region, India. 

Figure 3 is a flow sheet describing the process of the invention. 
Detailed description of the invention 

The present invention relates to solving the governing stochastic heat conduction 
equation to obtain the mean and variance in the heat flow fields. The governing heat 
conduction equation in steady state is expressed as 



The present invention relates to solving the governing stochastic heat conduction 
equation to obtain the mean and variance in the heat flow fields as shown in Figure 1 . 

This heat conduction equation is mostly solved under the assumption of a constant 
thermal conductivity or temperature dependent thermal conductivity. Both these assumptions 
require one to assume realistic values of the thermal properties in the crust and the upper 
mantle. However, due to the heterogeneous nature it becomes difficult to estimate realistic 
values of these parameters. One way to overcome this problem is to consider the thermal 
parameters to be a realization of a random process. 

In the present invention the radiogenic heat production is assumed to be exponentially 
decreasing function and the thermal conductivity is considered to be a random parameter 
which is expressed as a sum of a deterministic component and a random component The 
stochastic heat conduction equation then reduces to 





(2) 



where T is the temperature ( *C) 



Ao is the radiogenic heat production at the surface ( /i W I m* ) 
D is the characteristic depth (m), 

K(z) « K + K'(z) is the thermal conductivity (Wlm'C) 
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which is expressed as a sum of a deterministic component and a random component ( K is 
the mean or the deterministic value and K (2) is the random component). 
Boundary conditions employed in this invention are 
Constant surface temperature 
5 T(0) = T 0 (3) 
and a constant mean surface heat flow 

f<-=Q, (4) 

The random component of the thermal conductivity is assumed to have a Gaussian 
colored noise correlation structure represented by 
10 E(K(z))=0 (5) 

E(K(z,)K(z,)) = aie^-l (6) 
where o 2 K is the variance is thermal conductivity, p is the correlation decay parameter (or 
1/p is the correlation length scale) and zi and zi are the depths. 

Following the procedure given in Geophysical J International, V 138, 1999, by 
15 Srivastava and Singh, the solution to the temperature field is obtained using the series 
expansion method as 

- (Q ' ~f' P) f(«-QK c ft)dS (7) 

In this invention the closed form analytical expressions for mean heat flow and error 
20 bounds on heat flow have been obtained The expression for heat flow with random thermal 
conductivity is expressed as 

q(2) = K*L + K'(z)^L (8) 
dz dz 

Using the solution of T(z) ( equation 15) and differentiating it by using Leibnitz rule 
for differentiation under integrals and then taking expectation on both sides of equation (8) 
25 and applying the property of random fields ( equation 5 & 6) the expression for the mean 
heat flow has been obtained as 

q(2) = (1 - C* )(Q, - A 0 D( 1 - e- l/D ) ) + C\ Q, e" (9) 

where C K = a K / K is the coefficient of variability of thermal conductivity. If C K = 0 in 

the above equation then the solution reduces to the deterministic solution. 
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From the above equation it is clear that the mean heat flow is controlled by the 
coefficient of variability of thermal conductivity and to a small extent on the correlation 
length scale. 

The random component of the heat flow can be expressed as 
q(z) = q(z) - q(z) 

- l&jpifcM d4 + «(*)(Q,-A 0 D + A 0 De"> ) (10) 

The covariance at zi and z 2 is obtained by taking expectation of E{ q (z,) q'(z 2 ) } and 
setting zj = zj = z the closed form solution to the variance in heat flow has been obtained as 
a] = q 1 x Tl + q 2 x T2 + q3 x T3 + q4 x T 4 + qS x T5 - q6 x T6 + ?7 x Tl - 98 x Ti + q9 x T9 

(JJ) 

where 

q\=C 1 K Al(\-pD) 1 

q2 = q3 = C 7 K (Q t -A 0 D)A o p{pD-l) 

q4=C 2 K (Q,-A 0 D) 2 p 2 

qS=ql = C* {Q, -A 0 D + A 0 De-"°)A 0 (1 - pD) 
96 = 98 = Cl (Q, - A 0 D + A.De-'^XQ-A.D) p 
q9 = C 1 K (Q,-A i) D + A 0 De"">) 1 
The closed form solution to the integrals are 



1 f(l-e- z « /0 ) (e 



1 ((!-<? 



(p + l/£) J 



-Ix/D 



(p + l/D)^ 2/£> 



1 „ fe 



-i(X7*l/D) _ e -2«/D ^ 



72= - 
P 



( (\-e-" D ) ( e -<"*' /0 > -I) 
l/D + (p + l/D) 



1 

P 



f (l-e" /D ) frr" -e- /p ) ' 
l/D + (p-l/D) ; 



T3 = 



1 



(p-1/0) 



f (l-g'"°) (V" -1) 1 
I/Z> + p 
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f O-e" /J ) ( 



T4= ^| 2 + £^ 



TS _ (e-" D -e-") 
(p-l/D) 



7-6 = i T7 = T5 T8 = T6 T9 = 7 

I 

The novelty in the present invention over the prior art methods resides in that the 

S method uses randomness in thermal conductivity structure to quantify errors in heat flow and I 

leads to closed form solutions to the mean and its error bounds on the subsurface heat flow i 

The present invention describes a method of solving the stochastic heat conduction 

equation to obtain the closed form solution to the mean and variance in the heat flow fields. j! 

Using the thermal conductivity to be a random parameter the stochastic heat conduction j 

10 equation has been solved and first the solution to the temperature field is built using a series j 

hi 

expansion method. Next the expression for heat flow is obtained and subsequently the ? 
expressions for the mean and variance in heat flow have been derived. jjj 
The aim of this invention is to provide exact formula for quantifying the error bounds jjj 
on the subsurface heat flow. * I 

15 EXAMPLES 

The subsurface heat flow along with its error bounds up to Moho depth has been '■ 
computed using the expressions obtained for mean heat flow and its standard deviation for " 
two examples (1) Synthetic Earth Model and (2) Latur region of India, which shows the jj 
efficacy of the present invention. Most of the present investigations compute the subsurface 

20 heat flow in a deterministic way where the solution are obtained on the assumption that the 

input controlling thermal parameters are assumed to be known with certainty As the Earth | 
system is highly heterogeneous some amount of uncertainty is bound to exist in the input W 
controlling parameters Hence one needs to quantify the errors in the dynamical behavior of * 
the system. In this invention the error bounds on the subsurface heat flows have been 

25 quantified for different values of the controlling thermal parameters. 

Example 1: Controlling input thermal parameters for a Synthetic Earth Model 
Numerical values of the controlling input parameters for a realistic Earth model are. 
Radiogenic heat production (A) 3.0 ( ^W/m s ) 

Characteristic Depth (D) 1 2 (km) 

30 Moho Depth (L) 50 (km) 



i 
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Surface Temperature (T 0 ) 0 ( m C) 

Surface heat flow (Q,) 65 (mW/m 2 ) 

Random thermal conductivity: 

Mean thermal conductivity AT 3.5 (mW/m 2 ) 

Coefficient of variability C k 0.3 

Correlation length scale 1/ p 25 km 

Example 2: Controlling input thermal parameters for the Latur region, India. 

Numerical values of the controlling input parameters, Tectonophysics, V 306, 1999, 
Roy and Rao. 

Radiogenic heat production ( A) 2.6 ( /ufV/m' ) 

Characteristic Depth (D) 12 (km) 

Moho Depth (L). 37 (km) 

Surface Temperature (T 0 ) 30 CO 

Surface heat flow ( Q s ) 43 r mW/m 2 ) 

Random thermal conductivity: 

Mean thermal conductivity K 3,0 ( mW /m 2 ) 

Coefficient of variability Cr 0.2 
Correlation length scale 1/ p 17 km 

The parameters associated with the random component are the coefficient of 
variability in the thermal conductivity and the correlation length scale The coefficient of 
variability in the thermal conductivity can range from 0 to 1 0 and the correlation length scale 
should be less than L i e the total length of the model For this numerical study we have 
considered 20% error in the thermal conductivity and about half the length of the model as 
the correlation length scale. If we wish to increase the errors in the input thermal conductivity 
and see its behavior on the system output we could increase the coefficient of variability in 
the thermal conductivity. 

Using these controlling parameters the mean heat flow has been computed using 
equation (17) and the standard deviation is obtained from the expression for the variance 
(equation (19)) Figure 1 shows the plot for mean heat flow along with its error bounds for 
example I and Figure 2 is the plot for mean heat flow along with its error bounds for the 
Latur earthquake region. From both the figures we see that the error bounds on heat flow is 
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maximum at the surface and decreases with depth. These error bounds are very important and 
these exact formula have not been given so far. 

In the present invention the expressions for mean heat flow and error bounds on heat 
flow have been obtained. The heat conduction equation is solved using a stochastic approach 
incorporating Gaussian randomness in the thermal conductivity and the solution is built 
using a series expansion method and closed form analytical solutions of mean and variance 
for heat flow fields have been estimated. These expressions are seen to be highly dependent 
of the controlling thermal parameters such as coefficient of variability in thermal 
conductivity, on the correlation length scale, on the depth and on the other controlling input 
parameters Using these expressions one can estimate the mean heat flow at any given depth 
along with its error bounds In most investigations the. heat flow is estimated in a 
deterministic way or the error are quantified using the Monte Carlo simulation method. The 
advantage of this approach is that closed form expressions for mean and variance in the heat 
flow have been obtained and can be directly used for estimating the errors in the heat flow 
field due to error in the thermal conductivity. 

It must be understood that the above description is intended to be illustrative of the 
invention Modifications are possible without departing from the spirit of the invention 
Advantages: 

1. The advantage of this invention is that exact formulaes have been given to quantify 
the error bounds on the subsurface heat flow 

2. The exact formula obtained for the mean and standard deviation in heat flow have not 
been so far given and will benefit the Geoscientific community dealing with Geodynamical 
problems This study can be used in quantifying the heat flow at any given depth along with 
its error structure for a conductive Earth model. 
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